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Abstract 

A simple perturbation description unique for all signs of curva- 
ture, and based on the gauge-invariant formalisms is proposed to 
demonstrate that: (1) The density perturbations propagate in the 
flat radiation-dominated universe in exactly the same way as electro- 
magnetic or gravitational waves propagate in the epoch of the matter 
domination. (2) In the open universe, sounds are dispersed by cur- 
vature. The space curvature defines the minimal frequency uj c below 
which the propagation of perturbations is forbidden. 

Gaussian acoustic fields are considered and the curvature imprint 
in the perturbations spectrum is discussed. 

PACS numbers: 98.80 Hv 
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1 Introduction 



Discovering the wave nature of scalar perturbations in the early universe has 
a long history. Watchful reader of Harrison's classical paper U can guess the 
wave equations out of formulae given there (Section 5.5). Trigonometric or 
Bessel solutions together with (^-dependence characteristic for flat perturbed 
universes appear in both classical |2|, [J and gauge- invariant theories ||. In 
the case of flat universe the gauge-specific wave equations are explicitly given 
by Sachs and Wolfe (0 see theorem pp. 76-77). The comprehensive phonon 
description of perturbations in the flat radiation-filled universe, together with 
the attempt to quantize them, has been formulated by Lukash |5| and contin- 
ued in its quantum aspect by others . The wave character is confirmed || 
in the original Lifshitz-Khalatnikov || formalism. Acoustic motions of the 
baryon-electron system after recombination have been noticed by Yamamoto 
at al. [|K|]. Some parallels between scalar perturbation dynamics and gravi- 



tational waves can be found in TT1, 112 



Controversies, however arise over the gravitational instability criteria, the 
gauge problems and the role of the space curvature. 

(1) In the 7] — > limit one can formally construct the growing and decay- 
ing solutions. Since these solutions are typically considered as the large scale 
approximation (a; — > 0) the structure formation is expected in scales greater 
than the sound horizon. Consequently the Jeans criterion is understood as 
the dispersion relation dividing perturbations into two classes: acoustic waves 
and gravitationally bound structures. No dispersion relations like that can 
be inferred from the exact solutions @]-[E|> ||. 

(2) As long as the results depend on the coordinate system (the gauge- 
specific solutions ||, |6], |8|] differ one from another) their physical meaning 
is a subject of dispute. Acoustic field deserves complete gauge-invariant 
treatment. 

(3) The problem of acoustic field does not seem to be solved properly 
in open universes, where most authors traditionally employ the flat space 
Fourier analysis, instead of Fourier expansions in the Lobachevski space. 

In attempt to clear those points, we propose a simple perturbation de- 
scription, which is unique for all signs of curvature, and based on the gauge- 
invariant perturbation formalisms (Sakai ||, Bardeen Q, Kodama and Sasaki 
[ |T3|j , Lyth and Mukherjee [nj], Padmanabhan [JI5] , Brandenberger, Kahn and 



Press flql , Ellis, Bruni and Hwang [17|, Olson [lq|, see also [T 
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Section |^ contains a brief recipe of how to reduce equations obtained 
in these theories to a single, second order partial differential equation (J3|). 
Differences between the formalisms occur to be of no importance here, and 
we obtain exactly the same propagation equation for all of them. We show 
how to transform this equation to the wave equation in its normal form. 

We obtain a general, "profile-independent" solution for the flat universe 
(Section |3|), without appealing to the Fourier transform. We demonstrate 
that the gauge-invariant density perturbation propagate in radiation-domi- 
nated universe in the same way as electromagnetic or gravitational waves 
propagate in the epoch of the matter domination. Eventually, we expand 
perturbations into planar waves, in order to discuss some basic features of 
the spectrum and the spectrum transfer function. 

In section f| we describe the sound propagation in open universes. We 
analyse the dispersive role of the curvature. The space curvature prevents 
perturbations of frequencies smaller than some critical u c from propagating 
in space, and systematically reduces the group velocity for others, when uj 
goes down to u c . 

Section |^ is devoted to Gaussian acoustic fields. We derive the spec- 
trum transfer function in the form suitable to estimate the role of the space 
curvature in the microwave background. 



2 Scalar perturbations in the early universe 

In the universe filled with highly relativistic matter the energy momentum 
tensor is trace-free. The dynamics of the scale factor a{rj) expressed as a 
function of the conformal time r\ is governed by 

T^ = --^- ) (a"(r ] )+Ka(r ] )) = (1) 

and yields 




We treat the curvature index K as continuous quantity and keep K explicitly 
in both the equations and solution, as far as possible. Traditional formulae 
can be recovered by setting K — ±1 or by the limit procedure K — > 0. 
Normalization recalls the constant of motion M. = p(r])a 4 (ri). 
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The perturbation equation expressed in the orthogonal gauge?], |18|, |I9| 
and parameterized by the conformal timef] 77 takes the canonical form (free 
of first derivatives) 

ff IhA 1 

;X(V, x) - -A^Xiji, x) - - ®AX( V , x) = (3) 



drf ' 3a 2 (rj) ' 3 

where ®A denotes the Laplace-Beltrami operator acting on orthogonal hy- 
persurfaces. 

Equation (|3|) can be easily derived from the Raychaudhuri and the con- 
tinuity equations (see the procedure in [TJj] or |I5|| ). It also can be recovered 
from the Sakai equation || (formula (5.1), K — > X), the equation for den- 
sity perturbations in orthogonal gauge (Bardeen's || formula (4.9), p m — > X, 
Kodama and Sasaki chap. IV, formula (1.5), A — * X, Lyth and Mukher- 
jee formulae (16-17), 5 -> X, Padmanabhan pf Eq. (4.88), 6 -»• X), 



the equation for gauge invariant metric potentials (Brandenberger, Kahn 
and Press |16| formula (3.35), $H/pa 2 X), the equation for gauge invari- 
ant density gradients (Ellis, Bruni and Hwang |T7| formula (38), T> — > X) 
or Laplacians (Olson [18|] formulae (8-9), as well as its extension to open 
universes []B5| formula (22)) after transforming these equations to conformal 
time (if parameterized differently) and employing the Helmholtz equation to 
restore partial form of the perturbation equation. Suitable changes of the 
variable names as indicated above (original X) are necessary. 
We introduce a new perturbation variable X 

X(r ? ,x) = -i y |-(a(r ? )X(r ? ,x)). (4) 

While X(j], x) satisfies (§) and a(rj) is given by (@) the perturbation vari- 
able X obey the wave equation in its normal form 

^X(r / ,x)-l(%X(r / ,x) = 0. (5) 

The time derivatives of the gauge-invariant inhomogeneity measures also 
form gauge-invariant variables. In this sense variable X is equally good as X. 

1 oi the gauge-invariant differential measures of inhomogeneity []17j 

2 In the orthogonal gauge the conformal time is defined as the integral r\ — J-^dt, 

where time t means the orthogonal time — the time parameter constant on orthogonal 

hypersurfaces. 
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However, time derivatives may be difficult to observe at the last scattering 
surface, and hardly represent physically meaningful aspects of the cosmic 
structure. The equation (|5[) plays only an auxiliary role, nevertheless is very 
useful. Formally, it describes the wave (massless field) propagating in the 
static space-time of constant space-curvature. The specific case of positive 
curvature (Einstein static universe) has been considered |20| in the context 
of quantum field theory on curved background. In our case spaces of zero or 
negative curvature are of particular importance. We will discuss both cases 
individually. 



3 Sound waves on the flat background 

When the space curvature vanishes the equation (^) reads as 

^- 2 x(v, x) - lx{ v , x) - i ^AX( V , x) = (6) 



and is essentially the same as the propagation equation for gravitational |TT 



12|| or electromagnetidj waves in the dust-filled universe 

d 2 „, , 2 



X( V , x) - -X( V , x) - ^AX( V , x) = 0. (7) 



drf 7] 

The only differences are that gravitational and electromagnetic waves are 
expressed by the tensor h^ u and vector respectively, and they propagate 
with the speed of light (c = 1), while the solutions to equation (||), represent 
scalar waves travelling with the velocity v — 1/ \^3. 

Now, the Laplacian ®A operates in Euclidean space. Equation (H) when 
expressed in Cartesian coordinates {x} is solved by an arbitrary function 
X = X(n-x - vr)) p| with v = l/y/3 and n-n = 1. However, to keep 



the linear approximation valid we require that X(n-x — vr]) to be limited 
throughout the spaceQ. 

Knowing the general solutions of the equation (||) we can return to ob- 
servables X. We look for the general form of X(rj, x) among the solutions of 



3 See formula (5.2.6) in |2l] after substituting g = if. 
4 Vr/ 3e e IR: Vx - e < X(n x - vr)) < e. 
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the equation (f|) 

X( v , x) = - (^J\' X(n-x - vr)')dr)' + F(x)) . (8) 



On the strength of (|6]) -F(x) is harmonic ' 3 'AF(x) = and must be constant 
if limited throughout the space of constant curvature |23[| . With no loss 
of generality^, we put F(x) = 0. Eventually, the general, spatially limited 
solution to the equation (|) is expressed by the integral 



1 rV — . 

X(j],x) — - rj' X(n-x — vr)')dr)' (9) 

77 •/ 



of an arbitrary, but also spatially limited function X(n-x — W7). The solution 
describes a wave having the time-dependent profile and travelling with the 
constant velocity v = l/v3. 

This can be easily confirmed by the Fourier expansion analysis. Indeed, 
for any real function X(n-x — vrj) expressed as 

X( V ,x) = J(hMv,x) + A*u*(^,x))dk (10) 

with 

u k {r],x) = ^=e l{k - x -^' ) (11) 
y2u 



corresponds X(rj, x) 

X(rj, x) = j\ (A k u k {r], x) + A* k u* k (r), x))dk (12) 
expanded into modes u k (r], x) 

u k ( V , x) = /i^e** = JL f 1 + — ) e^'*"^ (13) 



The frequency obeys the dispersion relation uj 2 = k 2 /3, the Fourier co- 
efficient A k = — — A& is an arbitrary complex function of the wave number 



5 The freedom to choose this constant is not different from ambiguity in the indefinite 
integral in (|^). Appropriate integration constants are traditionally tuned to give the 
perturbation spatial average equal to zero. 
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k, and Uk is obtained from (|9|) after substituting X = u^. Modes Uk like u& 
form an orthonormal base in the function space with the Klein-Gordon scalar 
product P0 |. Both and Uk form travelling waves, but only u^ have their 



absolute value constant in time. Therefore, the generic perturbation X(rj, x) 
is composed of plane waves Uk of decaying amplitude, which perfectly agrees 
with Sachs and Wolfe results (|| pp. 76-77) obtained in alternative perturba- 
tion approach. Waves move with constant velocity v = 1/ v3 independently 
of their length-scale. Short-scale and long-scale perturbations do not form 
different classes of solutions. 

In the theory appealing to stochastic processes the initial perturbation 
is given at random at the end of the quantum epoch rji > 0, and develops 
gravitationally according to (|5[) in the interval 77 > 77^. Therefore, solution's 
singularity at 77 = is purely mathematical fact with no physical conse- 
quences. 



4 Sound waves in the curved space 

While decomposing perturbations into Fourier series in flat or opened uni- 
verses we should respect some specific effects caused by the curvature. This 
particularly refers to open universes (the Lobachevski space), where the or- 
thogonal expansions exist only for a class of perturbations with sufficiently 
short-scale autocorrelation [p4| . To expand the others one needs supple- 
mentary series (supercurvature modes [25, 26, 27, [28|]) of non-orthogonal^ 
solutions to the Helmholtz equation, which are numbered by imaginary wave 
numbers k G [—*,*]. 

Let us adopt spherical coordinates {r, 9, <fi} as more appropriate for curved 
maximally symmetric spaces. In the same manner like in the scalar field 
theories |2{| the density perturbation expands as 



X(r], r,6,(f)) = J2 J ( A kim Uki m (v, r, 8, <p) + A* klm u* klm (r), r, 9, <p))dk (14) 

lm 

where modes u klm (r), r, 9, 4>) = ^,K{ri) Y kim{r,0A) are expressed by hyper- 
spherical harmonics Y k i m (r,9,(j)) and time-dependent amplitude ^uj,k{.v) ful - 

6 In the sense of the Klein-Gordon scalar product. 
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filling the time-equation (obtained by separation from (|3|)): 

d 2 , , (k 2 -K 2K 
dr] 



sin (yKr/) 



We find solutions to (yj|) in the exact form as 



V2' 



u 2 -K 



l + VK 



cot(v Kr\) 



-•ILxJTj 



'luJ 



(15) 



(16) 



and their complex conjugates. Solutions ^uj,k(v) approach Huirj) (eq. [13]) in 
the K — > limit. The frequency to and the wave number are related to each 
other by the dispersion relation 



u(k) 



Vk 2 -K 



(17) 



which can be obtained by simple substitution of (0) into ( |15"D and perfectly 
agrees with the dispersion relation obtained for the variable Y on the strength 
of equation (||) (compare |20] chapter 5.2). 

Functions Y k i m (r, 9, <fi) solve the Helmholtz equation [ED 



(3) AY klm (r,9 



\k 2 - K)Y klm (r,#, 



(18) 



and can be split into the radial part 11^ and the two dimensional spherical 
functions Yi m (-d, ip) 



Y k i m (x,^,if) = U k i(x)Yi m ('d,ip) 
Solutions to radial equation 



dx 

are given by 



2 n w (x) + 2cothx^n w (x)- f A + ^2^ 



(19) 



U kl { X ) = (20) 



N kl 



N,,T1 



kl^kl 



.n=0 



-1/2 



l+l 



U kl = (k 2 smh X ) 1 



k sinh x d(&x) 



cos(kx) 
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The lowest multipole solutions 



n 



AO 



1 / sinfcx" 
sinhx y k 



~ if. smkx\ 

Ilfci = — r— - cos kx + coth x — : — 
smnx y k J 

= — ; — ( —3 coth x cos kx + (3 coth 2 % — A; 2 — 1) S " — ) 
sinhx V k J 

are enough to demonstrate properties of both two series of hyperspherical 
harmonics. For real wave numbers (k 2 > 0) the U^i (consequently also 
Ykim{r, 9, 0)) functions oscillate in space. They form an orthonormal basis in 
the sense of the scalar product (/1I/2) = / /i/2i/9^ x ' As proved by Gelfand 
and Naimark they are complete to expand square integrable functions in 



the Lobachevski space [30, [JT]. For imaginary wave numbers contained in 
the interval — 1 < k 2 < 0, the Hfci (and Yki m (r,6,(j))) functions build sup- 
plementary series. These functions are regular, limited but strictly positive 
throughout space, so they are not orthogonal]]. The supplementary series is 
redundant for expansion of square integrable functions. Nevertheless, this 
series is necessary to expand weakly homogeneous stochastic processes in the 
Lobachevski space |25], p6f . 

In this way in the open universe one obtains two types of hyperspherical 
harmonics Yki m (r, 9, 0) and consequently two types of modes Uki m (r], r, 9, </>). 
Modes Ukim{f] ', t ; 9 , 0) with real k are orthogonal by means of Klein-Gordon 



scalar product and expand waves of square integrable profile. Modes 
Ukim{f]i r i ®i 0) with — 1 < k 2 < form "waves of infinite length-scale". Both 
types of modes may contribute to the spectrum of randomly (or quantum) 



originated inhomogeneities p7| , 



7 Modcs with k = ±i are constant throughout space. One can subtract them by suitable 
changes in the background metric. Other modes with k G (— i, i), although positive every- 
where, decrease with distance strong enough to keep zero mean value. For instance, for 
the spherically symetric perturbation with the density excess given by p(x) — H±q(x) = 

2csch(v)sinh (%) the mass-to- volume ratio ~p(r) = — ^ Q — ; — ^ X tends to zero with 

\2J P\ > i7r J sin h 2 ( x )d X 

volume tending to infinity: p(r) — > ^J^h^r)^ — * 2 csch(r) sinh (|) — > 0. No redefini- 
tion of the background can absorb perturbations like that. 
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The density perturbations propagate in the open universe in different 
manner than the scalar fields or gravitational waves do. Acoustic waves 
of different length-scales propagate with different velocities. Indeed from 
relation ( |P7D we can infer both the phase and group velocity of sound in the 
form 

and 

* w = T^ k) = vmw (22) 

The group velocity decreases with the wave number k, to vanish completely 
at the k — > limit. The condition k = determines the critical frequency 
co>(0) = l/y/3, below which the wave propagation is forbidden. Therefore, the 
acoustic travelling waves are composed of the principal series modes. The 
supplementary series build 'global' standing waves of supercurvature scale. 



5 Gaussian acoustic field 

The generic acoustic field described by (||) is composed of waves travelling 
in different directions. Clearly, this property also refers to solutions of the 
equation The mechanism, which creates initial small perturbations is 
expected to be of probabilistic nature (thermodynamic or quantum fluctua- 
tions), therefore the evolution of linear structure is usually expressed in the 
language of stochastic processes. The homogeneity of this stochastic process 
reflects the universe homogeneity. Weakly homogeneous processes have their 
Fourier expansions 

X( V ,x) = |(A fe u fe (7 7 ,x) + AX(^x))dk (23) 

and consequently 

X{r], x) = J {A k u k {r], x) + A* k ut{r), x)) dk (24) 

where the coefficient A k are random variables of k and the integral has a 
stochastic sense [32]. In a generic Gaussian fieldf\the expectation values for 



8 Allen and collaborators |33[ ] distinguish between stationary random fields fulfilling 
(|25|-p6|) and the squeezed random fields, which violates (vw. Stationary processes have 
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Ak fulfill 

E[A k A* kl ] = V k 6 kV , (25) 
E[A k A k <] = 0. (26) 

Vk is defined as the field spectrum. The first relation guarantees the statis- 
tical independence of waves of different wave numbers, the second — says 
that no particular phase is preferred. Waves moving in different directions 
are statistically independent. 

The temperature fluctuations at the last scattering surface draw our at- 
tention to spatial correlations of Sp/p measured at the instant r] = r] r . In the 
flat universe the two-point spatial autocorrelation R(h) ]24| of the field X 
given by (|T^) can be expressed as: 

R{h) = ^-f E[X{x, r/)X[x + h, r])]5(h-h - 1) clh (27) 

2u k u* k V k exp(zk-h)5(h-h - 1) dk dh (28) 

rod 

/ A7ik 2 j (hk)2p^V k dk (29) 
Jo 

"f^^Wo* < 30 > 

where jo is the spherical Bessel function, and p k stands for the space spectrum 
of the density perturbation at a given moment rj. Following Peebles |34 



1 

47T 

roo 



WC 



define the transfer function T u {rj) = 2//^//* , which converts the time-invariant 
field spectrum V k into the space spectrum p k {rf). 

p k ( V ) = TMV k = 2p^V k . (31) 

The formula for the space spectrum p k (r]) splits into two factors: 1) T{rj) - 
describing the role of gravity, and 2) V k coming from other interactions and 
rendering their probabilistic nature prior or during radiational era. We do 
not discus any specific form of A k or V k in this paper. We assume, however, 
that A k enables one to construct small perturbations, and in particular does 
not cause divergences in Fourier integrals. Employing flT3"D one easily founds 




(32) 



their precise meaning in the framework of the stochastic theory not equivalent to ( |25[]2q ) . 
We prefer to name these random fields generic to stress analogy to generic classical fields. 
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As seen from ( p2|) the evolution of each mode depends on the product {ujrj) 2 . 
The contribution from modes much larger than the horizon scale urj <C 1 
strongly decreases with time, while the perturbation well inside the horizon 
ujr] ^> 1 keep constant amplitude. This property confirms the stability of 
Robertson- Walker symmetry against the generic (both classical and stochas- 
tic) large-scale density perturbations. 

In the same manner one can express random acoustic fields in the open 
universe 

X(rj, r,9,(j)) = Y, J (A klm u k i m (r], r, 6, <p) + A* klm u* klm (r), r, 0, <p))dk. (33) 

lm 

According to Yaglom's theorem, for weakly homogeneous processes integra- 
tion runs over both series principal and supplementary. The autocorrelation 



function now reads 26 



R(h)= I Anep^- Pk ( V )dk (34) 

J KSIII ft 

B+U[0,t] 

and the transfer function determined from formula ([16]) and (|3T| ) takes the 
form 

TM = -(l+ 2 l - . f ) (35) 
lo \ to 2 — K smh {yKvj) J 

We rewrite T u (rj) as a function of the energy density 



to demonstrate that the curvature modifies substantially the space spectrum 
p k only in the low frequency limit (supercurvature modes). Therefore, one 
may expect to find the curvature signature mostly in low multipoles (dipole, 
quadruple) (^7|. To extract this geometrical effect the knowledge of the field 



spectrum V k is indispensable. 

Expansions (|i~2"|) are typically employed in the gravitational waves theory 
and in the scalar field theory |2^| , while the density perturbations theo- 
ries traditionally solve, basically the same, propagation equation in terms of 
the Bessel J3/2 and Neumann N 3 / 2 functions. The J3/2 and N 3 / 2 are identified 
with "growing" and "decaying" modes respectively, according to their limit 
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behaviour at r\ = 0. Since the transition from {u, u*} basis to {J3/2, N 3 / 2 } is 
the unitary transformation, both representations are equivalent. When the 
"decaying" mode is rejected 



(a "standard practice" in cosmology — see 
comments in |36|) the unitarity is broken down and solution space is trun- 
cated to the space of standing waves. Then the acoustic field is in a highly 
"squeezed state" and consequently the characteristic peaks in the transfer 
function appears [|35|. 



6 Remarks on scales and observables 

There is a substantial difference between the dispersion on curvature, we de- 
scribed above, and "the curvature imprint" in CMBR spectrum, anticipated 
by the acoustic peaks hypothesis. The last hypothesis claims that the early 



perturbed universe was overdominated by stationary waves po| . One may 
justify this assumption by appealing to squeezing phenomena at the transi- 
tion from deSitter phase to the radiation dominated epoch. (In a transitions 
like that the large-scale stationary gravitational waves are generated (see 



57]). The same refers to massless scalar field.) The dominance of standing 
waves with specifically correlated phases should exhibit a series of peaks in 
the CMBR spectrum. Positions of these peaks are sensitive to details of the 
universe dynamics (fi,A, etc.). 

On the contrary, the dispersion effect described in the section [| has 
strictly geometrical character. It comes directly from the wave equation 
in the radiational epoch and does not depend on universe's past (evolution 
prior to radiational era is irrelevant here). The radiation-filled universe be- 
come "opaque" to sound waves greater than the curvature radius, whatever is 
the sound origin. In particular, no additional mechanism preferring standing 
waves at the beginning of the radiational era is needed. On that account the 
dispersion might form a reliable curvature tracer, provided it is observable 
at all, i.e. 

1) the space scales of supercurvature perturbations must be "small enough" 
to fit well in the observable part of our universe, and 

2) observational data should be complete enough to distinguish between 
standing and travelling waves at the last scattering. 

Answer to the first question is relatively easy and it was already formu- 



lated in the literature in terms of multipole decomposition |27J. We repeat 
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the same result below by use of simple geometrical consideration. Let us 
assume we live in the open (Q = 0.2) universe which is presently dominated 
by matter (p = 0) . We see the last scattering surface at some r\ r with the 
redshift z r = 1000. For sake of simplicity (following |37]]) we assume an in- 
stant transition from radiational epoch (p = p/3) to the galactic era (p = 0), 
which occur just at the last scattering moment rj — rj r . In such a universe 
model the scale factor evolves as 




a{rf) = 

and the radius of the visible universe Xr can be easily expressed as a function 
of redshift z r and the cosmological parameter Q 



z r ) = 2 arccoth — — 2 arcsinh 




Setting typical values Q = 0.2 and z r = 1000 one obtains Xr = 2.76. The 
equator plane "draws" on the last scattering surface a circle of the perime- 
ter I = 27rsinh(x r ) = 49.5, consequently, the curvature radius (in this units 
equal to one) takes a = 360// = 7.28 degrees on the sky. Let us consider 
now a spherically symmetric density perturbation (on the last scattering sur- 
face) described by the k — hyperspherical function Yqoo = x csch(x). The 
Yoi m functions, with arbitrary I and m, form a boundary between the prin- 
cipal and supplementary series, and can be understood as the "shortest" 
supercurvature modes. Since these functions are positive everywhere, we 
express the perturbation length-scale as the half-magnitude width ly2. For 
spherically symmetric mode looo it is roughly /i/ 2 = 2.2. This corresponds 
to 32° on the sky. Equator intersects about 10 patches of that size. The 
supercurvature perturbations contribute mainly to the lowest multipoles in 
the CMBR spectrum (roughly / < 5), but the visible part of the universe 
is large enough to produce the curvature effects. Another problem is how 
to determine from the CMBR data, which waves are stationary and which 
of them are travelling ones. The oscillation time scale for waves close to 



the critical frequency tu c = yl/3 is of the order of the universe age, thus 
no direct observation of their temporal behaviour can be done. The same 
refers to gravitational squeezed waves (see H33|). On the other hand, the 
information we need is "hidden" on the last scattering surface, and the "only 
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task" is to read it properly. For standing perturbations the density and the 
velocity fields are strongly correlated |38|. Density extrema coincide with 



the expansion extrema throughout the entire spaceg. This means that even 
in random acoustic fields the density and the velocity perturbations loose 
their statistical independence in scales comparable with the space curva- 
ture and consequently would gradually correlate on the sky when the angle 
scales outgrow ten degrees. The key to solve this problem is to find a sec- 
ond independent observable, which would help us to separate the potential 
and Doppler contributions to the CMBR temperature fluctuations. We can 
hardly propose a definite candidate at the present stage of observations, but 
both the polarization measurements and the large scale flows analysis seems 
to be steps in the right direction. 



7 Summary 

The gauge-invariant analysis confirms that density perturbations in the ra- 
diation dominated universe form a field of acoustic waves. In the flat uni- 
verse the density perturbations of all length-scales move with the same sound 
speed v — 1/ V3. Short and long perturbations do not form different classes 
of solutions. Perturbations' velocity is independent of the wave number, and 
in particular is the same for subhorizon and superhorizon inhomogeneities. 
Although the gauge-invariant theory confirms the fundamental properties 
of acoustic field, which have been known from gauge-specific descriptions 
[|j i? H) the propagation equations are not the same. An outstanding prop- 
erty of the gauge-invariant description is that propagation equation for sound 
in radiational era are identical with the propagation equations for gravita- 
tional or electromagnetic waves in the matter dominated universe. 

In the open universe perturbations evolve in a more complex manner. 
The negative space curvature causes the dispersion of acoustic waves. The 
universe geometry determines the minimal frequency for traveling acoustic 
waves in the similar way as the geometry of the wave conductor determines 
the minimal frequency for waves propagating inside. The critical frequency is 
related solely to the space curvature (not to the Jeans length-scale). Below 
this frequency perturbations form standing waves of supercurvature scale. 

9 Considerably more complex coincidences may be expected when multifluid models are 
taken into account. These cases need independent investigations. 
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In the radiation dominated universe the distinction between travelling and 
standing acoustic waves strictly coincides with the division into sub curva- 
ture and supercurvature inhomogeneities. Supercurvature standing waves 
are generic solutions. 

As commonly expected, the spectrum transfer function depends on the 
universe geometry, but differences are essential only in the large scale limit. 
In the subcurvature regime generic Gaussian acoustic field evolve like the 
acoustic field in the flat universe. Significant curvature effects may appear 
in supercurvature scales — lowest multipoles in the MBR temperature map. 
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